Effects of fluctuations on the initial eccentricity from the Color Glass Condensate in 

heavy ion collisions 
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We introduce a modified form of the Kharzeev-Levin-Nardi (KLN) approach for nuclear collisions. 
The new ansatz for the unintegrated gluon distribution function preserves factorization, and the 
saturation scale is bound from below by that for a single nucleon. It also reproduces the correct 
scaling with the number of collisions at high transverse momentum. The corresponding Monte 
Carlo implementation allows us to account for fluctuations of the hard sources (nucleons) in the 
transverse plane. We compute various definitions of the eccentricity within the new approach, 
which are relevant for the interpretation of the elliptic flow. Our approach predicts breaking of the 
scaling of the eccentricity with the Glauber eccentricity at the level of about 30%. 

PACS numbers: 12.38.Mh,24.85.+p,25.75.Ld,25.75.-q 



I. INTRODUCTION 

The elliptic flow V2 is one of the most important ob- 
servables in high energy non-central nucleus- nucleus col- 
lisions at RHIC It is very sensitive to both the ini- 
tial condition in the overlap zone, and to the subsequent 
evolution of the hot matter created in these collisions ■ 
Predictions from ideal hydrodynamics with an initial con- 
dition from a Glauber type model and with a short ther- 
malization time t < 1 fm/c are compatible with the el- 
liptic flow data at low transverse momentum at RHIC Q 
energies. This is one of the strongest indications of early 
thermalization and nearly perfect-fluid evolution of the 
quark-gluon plasma at RHIC. 

It was found in Ref. [J that ideal hydrodynamics 
with the initial condition taken from the Color Glass 
Condensate (CGC) approach Q overestimates the mea- 
sured elliptic flow at RHIC. Since hadronic dissipative 
effects were taken into account by a hadronic transport 
model @ ; additional dissipation in the quark-gluon phase 
would be required. The reason for the strong flow is that 
the KLN approach Q predicts a larger eccentricity of the 
initial gluon distribution in the transverse plane than the 
Glauber model [1, [l]. A recent calculation within the 
classical Yang-Mills (CYM) approach on the lattice Q 
gave somewhat smaller values than predicted by the KLN 
model. Nevertheless, the eccentricity of the "overlap" re- 
gion in coordinate space obtained from CYM still exceeds 
that predicted by the Glauber "wounded nucleon" model 
for soft processes. 

The eccentricity e is defined in the transverse plane 
i'± = ii"x,i"y) perpendicular to the beam axis by 



(1) 



The average may be taken with respect to either energy 
or particle density (or entropy, if rapid thermalization is 
assumed) . In the ideal hydrodynamical limit fl^ , elliptic 



flow only depends on the density and on e. Therefore, it is 
very important to know the initial conditions to extract 
properties of hot and dense matter created in nucleus- 
nucleus collisions. 

This paper is a continuation of Ref. , where we com- 
puted the eccentricity in the initial stage of heavy ion 
collisions within the KLN approach, and showed that 
the large eccentricity predicted by the fcj^-factorization 
approach is generic, as it does not depend on the details 
of the unintegrated gluon distribution function (uGDF) 
such as the presence of leading-twist shadowing and an 
extended geometric scaling window [ll|, [13, • 

The predictions of the CGC framework are consistent 
with RHIC data for the centrality, energy and rapidity 
dependence of the multiplicity of charged hadrons if the 
square of the saturation scale (averaged over the trans- 
verse plane) is defined to be proportional to the number 
of participants 0, H, [13, In this paper, we shall 

consider in more detail the relation between the differ- 
ent definitions of the saturation scales (see also the short 
note in Ref. [3). 

The eccentricity could of course fluctuate for many 
reasons, and it was recently argued that event-by-event 
flow fluctuations can affect the elli ptic flow obtained via 
correlation methods [3, [13, [H, 1X9^ . In order to ex- 
plore such fluctuation effects in the CGC picture, we em- 
ploy a Monte-Carlo implementation of the KLN approach 
which accounts for fluctuations of the high-rapidity "hard 
sources" (nucleons) in the transverse plane. This leads 
naturally to a factorized version of the KLN approach, 
where the uGDFs depend only on the properties of a sin- 
gle nucleus. We analyze the effects of the fluctuations on 
the eccentricity by using different definitions of the ec- 
centricities [IB, [13, [H, [l3| • This is relevant to the study 
of the effects of fluctuations on V2 [13, [Ml, [13, [11 • 
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II. GLUON PRODUCTION IN THE KLN 
APPROACH 

In the fcj^-factorization approach [l^l , the number dis- 
tribution of produced gluons is given by 



dNr, 



X {p± + k±f/A) (j)B{x2, (p_L-fc±)V4) 

with Nc — 3 the number of colors. Here, p± and y denote 
the transverse momentum and the rapidity of the pro- 
duced gluons, respectively. The light-cone momentum 
fractions of the colliding gluon ladders are then given by 
2^1,2 = P-L exp(±j/)/-yi, where ^/s denotes the center of 
mass energy and y is the rapidity of the produced gluon. 
We set p™'^^ such that the minimal saturation scale 
(58,111111(3^1,2) in the above integration is Aqcd = 0.2 GeV 
(see, however, below). The KLN approach 0] employs 
the following uGDF: 



1 



Ql 



as{Ql) max(g2,A:2) 



(3) 



where Qs denotes the saturation momentum at the given 
momentum fraction x and transverse position r±. The 
overall normalization is determined by the multiplicity at 
midrapidity for the most central collisions. 

In Refs. d, H, 0, Q , the saturation scale for nucleus A 
is taken to be proportional to the density of participants, 
n^3^j.j(rj^), which in turn is obtained from the thickness 
functions Ta and Tb'- 



(1 



(1 



b/2) 



b/2)/B)^) .(4) 



b = (6, 0) is the impact parameter vector in the transverse 
plane and aNN the nucleon-nucleon inelastic cross sec- 
tion. 

The density of participants is not a universal quantity 
which depends only on the properties of a single nucleus. 
In the following we propose a modified form of the KLN 
approach where factorization is manifest and show that 
the Ql ~ npart definition is in fact a good approximation 
for computing the multiplicity of gluons. Other quanti- 
ties such as the energy density, however, do behave some- 
what differently. 



HI. FACTORIZED KLN FORMULA FOR 
NUCLEUS-NUCLEUS COLLISIONS 

Generically, it is reasonable to define the squared sat- 
uration scale to be proportional to the density of nucle- 
ons. However, a problem arises near the edge of a nucleus 
where the density is small. There, the thickness function 
Ta is small: only for some configurations of nucleons do 
we actually find a nucleon at that position. Let us de- 
note this probability of finding at least one nucleon at 




FIG. 1: (Color online) Schematic sketch of different construc- 
tions of the uGDF. The full line corresponds to the uGDF of 
a single nucleon, the other lines to that of a nucleus at its sur- 
face. The dotted line corresponds to the the definition from 
Eq.® while the dashed line is for Eq. ©. 



a given transverse coordinate as pA- For configurations 
where indeed there is one single nucleon at the given , 
the uGDF should be that of a nucleon: (f> = (t){Q'\^). 
However, the average 0^(rj^) would then be 



rather than 



=PA(l){Qln) 



cj)A = (j){Ql^nPA) 



(5) 



(6) 



This is illustrated in Fig. [T] The essential point is that 
one has to average the uGDF itself, and not its argument 
Ta- Expressing the uGDF in terms of Ta, the relevant 
variable would be Tia, which is the density of nucleus A 
averaged only over those configurations with at least one 
nucleon at a given transverse position. This quantity can 
be obtained by averaging over the thickness functions of 
i> 1 nucleons, 



TiA = {t 



A}i>l 



E.=iPA(i) tAii) _ Ta 



Y.i=iPA{i) 



PA 



(7) 



where i^(i) is the thickness function for i nucleons in a 
row. This allows us to construct an uGDF in the follow- 
ing way. 



PA 



Ta 

PA 



(8) 



Here, and in the following, we use a short-hand nota- 
tion for the dependence of cj) on some density p: 4>{p) = 
(f){Qs{p))- The uGDF ^ respects factorization, since it 
depends only on the properties of a single nucleus. The 
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(9) 



saturation scale is then parameterized as 

n2 ( ^ nr. .r2 ( TA{r ^) / pA{r ^ )\ (Om 
Qs.Ala;,^-^) = 2GeV I — 

PA can be taken from the Glauber definition of ripart^ 

PA - (1 - (1 - STa/A)^). (10) 

In this case, S* is a sampling area and not the nucleon- 
nucleon inelastic cross section as in the Glauber model. 
Nevertheless, these quantities are related to each other, 
since 5* should be on the order of the area of the nucleon. 
One could fix its value in such a way that one recovers 
the saturation scale of a nucleon in the limit Ta — *■ 0, 
since 



lim 

Ta^O pa 



1 

5^' 



(11) 



For our simulations we will however use the value of the 
inelastic cross section at RHIC energy, since then the 
results are more directly comparable to the original KLN 
ansatz. The saturation scale of a nucleon is then equal 
to Ql„ = 2.0 GeVV(4.2 x 1.53) 0.31 GeV^ We refer 
to this approach as the factorized KLN (fKLN) model. 

Next, we explore the relation to the original KLN 
ansatz, where the saturation momentum is defined via 
n-paxt- The transverse density of produced gluons can 
roughly be ex]3ressed analytically in terms of the satura- 
tion scales as 



dy<Pr\ 




(12) 



where, Q^^max and Qs.m\n denote the larger and the 
smaller value of the two saturation scales in opposite nu- 
clei at any fixed position in the transverse plane. It is 
easy to see that this function is homogeneous of order one 
in both ^ and ^, or in the corresponding densities 
which appear in the definition of the saturation scales. 
The fcj^-factorization formula with the new definition of 
the uGDF (fKLN approach) now takes the form: 



dN„ 



rpi_ 
Pi 



d^k±_pAPB(t>{TA/PA)4'{TB/Pl 



(13) 

Due to the homogeneity, we can move the prefactors 
in eq. ([T3| into the arguments of the uGDF. If S is 
the nucleon- nucleon inelastic cross section ctjvat, we have 
J^part.A = TaPb- This leads to the original KLN defini- 
tion of the uGDF, where the saturation scales depend on 
the participant densities: 



dN„ 



dyd^ri^ 



fp± 

Pi 

d'^pi_ 

Pi 



d^k^4>{pAPBTA Ipa ) (j>aPb Tb/pb) 



fc_L (?!)(npart ,A ) </>('^part ,B ) 



(14) 
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FIG. 2: (Color online) Centrality dependence of the multiplic- 
ity for the different assumptions of the uGDF. Eq. (|6} gives 
a significant overestimate of the multiplicity in the peripheral 
collisions. The original KLN ansatz, where Ql ~ Wpart is very 
similar to the fKLN result which uses Eq. (|13p . Data from 
PHOBOS fH. 
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FIG. 3: (Color online) The initial eccentricity e at midrapidity 
as a function of the number of participants for 200 A GeV 
Au+Au collisions: KLN/e (KLN/n) corresponds to the KLN 
result using the energy (number) density as the weight in the 
definition of the eccentricity. fKLN results follow the same 
notation. The KLN results without a saturation scale cutoff 
are shown by symbols. 



The expression (jl3p involving the new definition of the 
uGDF can now be integrated numerically. The centrality 
dependence of the multiplicity at midrapidity is shown in 
Fig. [21 assuming that the hadron yield is proportional to 
the yield of produced gluons. For comparison, we also 
show the result corresponding to the definition ^ Ta 
in the fcj_-factorization formula, which clearly overshoots 
the data for peripheral collisions. On the other hand, 
the fKLN approach is rather similar to the original KLN 
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approach. This is due to the fact that the single-inclusive 
gluon cross section is homogeneous of order one in both 
arguments, as discussed above. 

In Fig. [31 we compare the eccentricities obtained from 
fKLN (eq. ^) and KLN. We plot two cases, where we 
use the transverse energy density dE±/dy(Pr± and the 
number density as a weight in the definition of the eccen- 
tricity. The eccentricity from fKLN is somewhat lower 
than the one from KLN, especially for more peripheral 
collisions. This is due to the fact that the transverse en- 
ergy distribution from the KLN approach is homogeneous 
of order 3/2, because of another factor of Qs,max US]- 



d'^r±dy 



Qs,minQs,max log 



Q 



2 

s,7nax 



)2 . 

■> s.7mn 



(15) 



As a consequence, the energy density distributions in the 
transverse plane from KLN and fKLN, respectively, are 
somewhat different. 

If one uses the number density of produced giuons 
rather than their energy density as a weight in the defini- 
tion of the eccentricity, one would expect that the fKLN 
form reproduces the KLN result, since the same argu- 
ments hold as for the multiplicity. However, we have 
found that the fKLN prediction for the eccentricity with 
the number density as a weight is almost identical to 
the above result (which employed energy density as a 
weight). This can be traced back to the above-mentioned 
cut-off Qs,min = Aqcd- In the fKLN model no such cut- 
off is needed as the minimal saturation momentum never 
drops below the saturation momentum of a single nu- 
cleon, by construction. Rather, the probabilities pA-, Pb 
of encountering at least one nucleon vanish at the sur- 
face of the overlap region. Even though it does not make 
much sense physically, technically we can set the mini- 
mal saturation scale in the KLN model to zero. Then the 
eccentricity from the KLN approach agrees well with the 
one from the factorizcd fKLN form as seen in Fig. [31 



IV. MONTE-CARLO IMPLEMENTATION OF 
THE KLN APPROACH 

According to the ideas of McLerran and Venugopalan, 
the gluon field of a sufficiently large nucleus can be viewed 
as a classical Yang-Mills field [20]. Particle production 
in high-energy hadronic collisions can then be calculated 
by solving the Yang-Mills equations in the forward light- 
cone with a boundary condition corresponding to clas- 
sical sources on its two branches. Real time solutions 
of the classical Yang-Mills fields on the lattice were ob- 
tained, for example, in Refs. [2^,[2^,[3^. Those numerical 
solutions arc usually averaged over many different initial 
conditions (configurations of the color-charge sources). 

Here, we calculate the thickness function used for 
determining the saturation scale with a Monte Carlo 
method. There are two advantages of doing this. First, 
the Monte-Carlo implementation naturally leads to a def- 



inition of the saturation momentum in terms of the thick- 
ness functions (which preserves factorization). Second, 
we can study the effects of event-by-event fluctuations on 
the eccentricity within the KLN approach. We only con- 
sider fluctuations of the positions of the nucleons, which 
were previously studied in Refs. [l^, [l3, [13, E| within 
the Glauber approach. We note that in Ref. [30|, a sim- 
ilar method was employed to impose a color neutrality 
condition for a each nucleon in the classical McLerran- 
Vcnugopalan model in order to simulate collisions of finite 
nuclei. 



A. Model 

Let us now construct a Monte Carlo model for gluon 
production for each configuration of nucleons within the 
colliding nuclei (denoted as MC-KLN). We first deter- 
mine the positions of A nucleons according to the Woods- 
Saxon distribution. Correlations among the nucleons are 
neglected at present but could in principle be taken into 
account. They tend to reduce fluctuations somewhat [3]| . 

The local density of nucleons, tA, is then obtaine d by 
counting nucleons within a tube of radius Tmax = ■^/S'/tt 
around a given transverse coordinate rj^: 



number of nucleons 



(16) 



The area S is the sampling area over which we count 
nucleons. The average (t^(rx)) reproduces the convo- 
luted thickness function TA{r±) with a sampling proflle 
/(^J_) = ^Q{r„iax - \r±\): 

{tA{r±)) = fA{ri_) = /rfVl r^(rl) /(r^-rl) (17) 



In the limit S' ^ 0, Ta equals Ta- This limit can be taken 
in a Monte Carlo by the oversampling ("test-particle") 
technique, whereby the sampling area S is reduced by 
some factor n and the number of nucleons is increased 
by the same factor; each nucleon then carries a weight 
1/n. (We found that oversampling by a factor of n 100 
gives stable results.) It is not appropriate to simply de- 
crease S. While this would reproduce the mean value of 
Ta, it would also change the fluctuations in the number 
of nucleons at a given position, since each nucleon would 
contribute n/ S to the density. With the oversampling 
technique each nucleon has a weight 1/n and therefore 
contributes \/S to the density, which is independent of 
n. However, correlations between neighboring points in 
the transverse plane are destroyed. Therefore, event-by- 
event observables cannot be obtained with this method; 
we employed it only to test convergence of our MC im- 
plementation to the "mean-field" results presented above 
and for the computation of averaged observables. In the 
following simulations, we choose S = <J\nc\ as in the mean- 
field fKLN approach above. 
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We then follow the same procedure to compute tB{r±), 
thereby generating two configurations for nucleus A and 
B, respectively, at a specified impact parameter. The 
saturation scale at a given transverse coordinate is then 
given by 



ix,r±) = 2GeV^ 



tAir±] 
1.53 



0.01 



(18) 



and similarly for nucleus B. It is clear that this choice 
of the saturation scale only depends on the properties of 
one nucleus, and therefore respects factorization. 

For each generated configuration we apply the k±- 
factorization formula ^ to obtain the distribution of 
produced gluons in the transverse plane. We then av- 
erage over many collisions to compute observables and 
their fluctuations. 

Finally, we establish the relation between MC-KLN 
and fKLN. The average uGDF which is effectively em- 
ployed in the MC-KLN model may be written as 

A A 

<I^A = 5Ip^(*)<^(*^(^)) = , (19) 

4=0 1=1 

where we used ^^(0) = 0. The fKLN uGDF may be 
obtained from the approximation where tA{i) is replaced 
by the average, tA{i) Tia- 

A 

(pA « ^PA(*)0(riA) =pa</'(Ta/pa), (20) 

i=l 

where J2i=iPA{'>') = PA was used. 



B. Results 
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FIG. 4: (Color online) Centrality dependence of the multiplic- 
ity in Au-I-Au collisions for ^/Inn = 200 GeV for the fKLN 
and MC-KLN models. All curves have the same normaliza- 
tion, we adjusted the result for KLN to the two most central 
bins. The data is from PHOBOS [H. 
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First, we show the centrality dependence of the multi- 
plicity from the MC-KLN model in Fig. |4] together with 
the result of the original KLN approach. For clarity, we 
have not adjusted the overall normalization to the most 
central events. The multiplicity as a function of the cen- 
trality from MC-KLN is very similar to the original KLN 
and to the fKLN approach. The overall normalization is 
about 5% lower. The reason is that the multiplicity as 
a function of the density has negative curvature (second 
derivative), so fluctuations reduce the multiplicity. 

In Fig. [5l we compare the eccentricity among the var- 
ious approaches. Since the fKLN result was obtained 
using a delta function sampling profile (by using Ta in- 
stead of Ta), we employed the oversampling technique 
with a factor n = 100 for the MC-KLN model. This way, 
the two models use the same sampling profile and the 
results are comparable. It can be seen that the eccentric- 
ity from MC-KLN is almost identical to the fKLN scal- 
ing, which is slightly higher than ricoii scaling. Note that 
any additional fluctuation (e.g. in the evolution of the 
BFKL-ladders) would not influence the standard eccen- 
tricity H]) , which is due to its deflnition; the denominator 



FIG. 5: (Color online) The eccentricity e in Au-fAu collisions 
predicted by the various models. fKLN and MC-KLN give 
almost identical results. The MC-KLN result was obtained 
using the oversampling technique with n = 100, which corre- 
sponds roughly to a delta function sampling profile, just as 
used in the fKLN model. 



and the numerator are averaged separately and both are 
linear functions of the weight. 

The eccentricity from the CYM approach 0] behaves 
similar to the number of collision scaling in the Glauber 
model. Note that the CYM approach in Ref. does 
not include effects of nucleon fluctuations. Hence, we 
conclude that the main reason for the reduced eccentric- 
ity as compared to the original KLN approach is the use 
of the thickness functions in the definition of the satura- 
tion scale as argued in Ref. Both, the CYM approach 
from Ref. @ and our fKLN approach, however, predicted 
larger e than the Glauber model for soft collisions where 
dN/dy(Pr± scales with the density of participants npart- 
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FIG. 6: (Color online) Various eccentricities for the MC-KLN 
model. No oversampling was used, which corresponds to a 
finite range sampling profile. 
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FIG. 7: (Color online) The standard eccentricity e and e-part 
of MC-KLN relative to the Glauber model result. 



C. Eccentricity Fluctuations 

We now turn to a comparison of various definitions of 
the eccentricity which are relevant for the scaHng prop- 
erties of the experimental data for V2 ■ Since we now con- 
sider cvcnt-by-event variables, we cannot use oversam- 
pling any more. This changes the results somehow, since 
the interaction range of the nucleon-nucleon interaction 
plays some role. In Appendix \^ we summarize how to 
include the interaction range in the numerical computa- 
tions of the number of participants and collisions. 

First we define the participant eccentricity Spart ac- 
cording to Ref. in order to account for fluctuations 
in the directions of the major axes of the overlap region 
and of the "center of gravity" , 

Spart = — 5 ' (21) 

where cr^ and Uy arc the RMS widths of the participant 
nucleon distribution projected on the r^, and ry axes, 
respectively, and Oxy = {xy) — (x) (y) . We also compute 

^{2} ^ V^(4lj (22) 

which is closely related to the event-plane method of the 
measurements of elliptic flow. For four-particle cumu- 
lants of azimuthal correlations [Toj we have: 

e{4} ^ {2{sl^,r - {el^,)y^\ (23) 

In Fig. [HI we show £part, £{2}, e{4} together with the 
standard e from the MC-KLN model, for Au-|-Au and 
Cu-I-Cu coUisions at = 200 GeV, using again the 
energy density as weight in the eccentricities. Compared 
to Fig. [51 the standard eccentricity is reduced somewhat. 



since we did not use the oversampling technique in this 
plot (see also Fig. [9] of Appendix A). As observed before 
H, [13, [ill in the Glauber model approach, e{2} deviates 
most from the standard definition, and the deviation is 
much larger for smaller systems. Since the source of fluc- 
tuations in this work is the position of the nucleons, the 
influence of the fluctuations on the results is very similar 
to that of the Glauber model approach. 

A detailed comparison to the Glauber eccentricity can 
be seen in Fig. [71 where we show the ratio of the eccentric- 
ity obtained by the MC-KLN model to the eccentricity of 
the Glauber model. In Au-|-Au collisions, we see an in- 
crease of the ratio eMC-KLN/£Giauber by more than 50% 
for central and 10% for peripheral collisions. Comparing 
to Cu-f Cu collisions, we notice a discrepancy of the ratios 
for a system size above 50 participants. Using the Epart 
definition, which accounts for fluctuations of the major 
axes, we observe that both systems follow the same curve 
and the enhancement is not more than 30%. Therefore, 
precise experimental data would be able to distinguish 
MC-KLN from the usual Glauber scaling by examining 
the scaling properties of the elliptic flow. 

So far we assumed that additional fluctuations (e.g. 
from the evolution of the BFKL-ladders) in the gluon 
production are negligible. While this is true for the stan- 
dard definition of the eccentricity, it does not necessarily 
hold for evcnt-by-evcnt definitions such as Epart, £{2} and 
e{4}. In these cases, the eccentricity might increase. We 
tested this assumption by including simple Poissonian 
fluctuations in the gluon number but found that £{2} 
increased by only about 0.02. 

V. DISCUSSION AND CONCLUSIONS 

We developed a simple model which includes fluc- 
tuations of the positions of the nucleons in the k±- 
factorization formula. It naturally leads to a gluon pro- 
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FIG. 8: The eccentricity as a function of the Qs,min low- 
density cutoff. 



duction cross section which respects factorization. The 
multiphcity as a function of centrahty does not change 
much as compared to the original KLN approach but the 
predicted eccentricity is somewhat lower. 

It should be noted, however, that the gluon distribu- 
tion calculated here corresponds to very early times on 
the order of ~ 1/Qs- Therefore, it is not necessarily iden- 
tifiable with the initial condition for hydrodynamics since 
it could be modified during the thcrmalization stage. A 
detailed analysis requires simulations of the time evolu- 
tion from the early production time l/Qs until thcrmal- 
ization. This could be done within parton cascade mod- 
els [H, 1131 which should be coupled to the Yang-Mills 
equations describing the soft gluons [U- This is beyond 
the purpose of the present paper. 

Our qualitative expectation is that the low-density sur- 
faces of the overlap zone may not thermalize (and should 
therefore not be included in the initial condition for hy- 
drodynamics), while the gluons produced closer to the 
center do. We can assess the consequences of such a sim- 
ple picture by increasing the minimal saturation scale 
Qs,'min{xi) which acts as a cutoff. This amounts to re- 
moving the very high-pj^ gluon jets as well as the low- 
density part of the bulk near the surface. The effect on 
the eccentricity is shown in Fig.[8l One observes that the 
eccentricity increases significantly with the cutoff. The 
fact that £ is flat up to Qs,min = 0.5 GeV is due to 
the construction of the fKLN uGDF: the saturation scale 
never drops below that for a single nucleon and smaller 
cutoffs therefore do not matter. Since no thcrmaliza- 
tion is expected to occur in p -\- p oi p + A collisions, 
the relevant cutoff for the hydrodynamic initial condition 
might be well above ~ 0.5 GeV. It would be interesting 
to study the hydrodynamical evolution of various initial 
conditions in order to check how sensitive the elliptic flow 
responds. 

Another advantage of the fKLN uGDF is that it re- 
produces the correct scaling of the multiplicity with 



"■ncoU = ctTaTb at high pt in the DGLAP regime. This 
is not the case for the original KLN ansatz. There, the 
high-p4 multiplicity also scales with Q^^ aQ^s ^ut with 
Qs,2 ~ "-part, one docs not recover scaling with TaTb- As 
we have shown, the naive modification to ^ Ta does 
not help, since the saturation momentum at the edge of 
a nucleus drops below that for a single nucleon, and the 
/cj^-factorization formula docs not reproduce the ccntral- 
ity dependence of the multiplicity. Our fKLN (and MC- 
KLN) uGDF solves these problems simultaneously: while 
it fits correctly the measured centrality dependence of the 
multiplicity, it does not break factorization, the nuclear 
saturation scale is bound from below by the saturation 
scale of a single nucleon, and it reproduces the correct 
ricou scaling at high pt- 

APPENDIX A: NUMBER OF PARTICIPANTS 
AND INTERACTION PROFILES 

1. Delta function interaction profile 

The densities of participants in the transverse plane 
are given by 

npart,A = (l - (1 - aTe/Bf^ (Al) 

npart.B = Tb (l - (1 - aTA/Af) (A2) 
"coll = (jTaTb . (A3) 

These definitions assume a delta function interaction pro- 
file /(r) = (7(5(r). In a Monte Carlo model, one can ap- 
proximate the delta function via the oversampling tech- 
nique: a a/n and A —f An, attributing a weight 1/n 
to each nucleon. 

2. Finite range interaction profile 

Next, we discuss how to calculate the number of par- 
ticipants and the number of collisions respecting finite 
range interaction. The thickness function is convoluted 
with the interaction profile, 

fAiri_) = J d^r'^l^^^^-^TAir'^) . 
Therefore, we have 

npart.A = Ta (l - (1 - ^tTs/S)^) (A4) 

npart^B = Ts (l - (1 - (JTa/A)^^ . (A5) 

The number of collisions can be derived in a similar way. 
Special attention is needed when specifying the location 
of the collision: ricoii = '^TaTb locates the collision at the 
participating nucleon from nucleus A and ricoii = gTbTa 
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FIG. 9: (Color online) Eccentricity for the Glauber model 
with number of participants scaling of dN/d?v±_ obtained 
with the following methods: delta function interaction profile 
(point nucleons), step function (hard sphere nucleons) and 
step function with sampling on a grid. 
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FIG. 10: (Color online) Same as Fig.|9]for number of collisions 
scaling. 



Glauber Monte Carlo, where the positions of the nucleons 
are sampled according to the Woods-Saxon distribution, 
and where the profile functions determine whether an 
inelastic interaction for a given pair of nucleons occurs. 
3. Monte Carlo Method on a grid 

Sampling the number of participants within the inter- 
action range corresponds to using the convoluted thick- 
ness functions for both nuclei. The corresponding defini- 
tions for the participant and collision densities are: 



ripart,A = Ta[1-{1- aTe/Bfj (A7) 
npart,B = Tb (l - (1 - aTA/A)^'^ (A8) 
"coll = ctTaTb (A9) 

Here the cross section is used in two ways: once for de- 
termining the interaction, and once for sampling the nu- 
cleons over the transverse area given by the cross section. 
Strictly speaking, this implies double-counting but this 
method provides good approximations to eqs. (jA4IIA6[) . as 
shown below. 

Figs . [9l and [TOl show the eccentricity for number of par- 
ticipants and number of collisions scaling obtained with 
different interaction profiles and sampling methods. The 
normal definition of Ta corresponds to a delta function 
interaction profile. It is reproduced by a Monte Carlo 
oversampling technique. A different in terac tion profile 
(in this case, a 0-fimction: /(r) = O^y/a/TT — r)) leads 
to formulas (|A4IIA6p . The corresponding Monte Carlo 
method (denoted as plain MC) is straightforward: the 
profile function determines the probability for an inter- 
action for each nuclcon-pair. The third method, Monte 
Carlo on a grid, samples the nucleon densities within 
the interaction range. It is very similar to plain Monte 
Carlo. However, the number of collision scaling is some- 
what smaller. 



locates the collision at the nucleon from nucleus B. The 
following definition 

ricou = I (TaTb + TaTb) (A6) 

uses the center-of-mass position and is symmetric under 
the exchange A ^ B. 

The definitions (|A4IIA6P agree with a straightforward 



Acknowledgments 

We are indebted to Adrian Dumitru for illuminating 
discussions and encouragement. We thank M. Gyulassy 
for useful comments, and T. Lappi and R. Venugopalan 
for interesting discussions. YN acknowledges support 
from DFG. HJD acknowledges support from BMBF grant 
05 CU5RI1/3. 



[1] B. B. Back et al. [PHOBOS Collaboration], Nucl. Phys. 
A 757, 28 (2005); J. Adams et al. [STAR Collaboration], 
Nucl. Phys. A 757, 102 (2005); K. Adcox et al. [PHENIX 
Collaboration], Nucl. Phys. A 757, 184 (2005). 

[2] J. Y. QUitrauh, Phys. Rev. D 46, 229 (1992). 



[3] for recent reviews see P. Huovinen, 'nucl-th/0305064'; 
P. F. Kolb and U. W. Heinz, nucl-th/0305084 
U. W. Heinz, nucl-th/0512051 ; T. Hirano, Acta Phys. 
Polo n. B 36, 187 (200 5); P. Huovinen and P. V. Ruuska- 
nen, |nucl-th/0605008| 



9 



[4] 

[5] 
[6] 

[7] 



[8] 

[9] 
[10] 

[11] 
[12] 

[13] 
[14] 

[15] 

[16] 
[17] 

[18] 
[19] 
[20] 

[21] 



[22] 



T. Hirano, U. W. Heinz, D. Kharzeev, R. Lacey and 
Y. Nara, Phys. Lett. B 636, 299 (2006). 
T. Hirano and Y. Nara, Nucl. Phys. A 743, 305 (2004). 
Y. Nara, N. Otuka, A. Ohnishi, K. Niita and S. Chiba, 
Phys. Rev. C 61, 024901 (2000). 

D. Kharzeev and M. Nardi, Phys. Lett. B 507, 121 
(2001); D. Kharzeev, E. Levin and M. Nardi, Nucl. Phys. 
A 730, 448 (2004) [Erratum-ibid. A 743, 329 (2004)]; 
Nucl. Phys. A 747, 609 (2005). 

A. Adil, H. J. Drescher, A. Dumitru, A. Hayashigaki and 

Y. Nara, Phys. Rev. C 74, 044905 (2006). 

T. Lappi and R. Venugopalan, nucl-th/0609021 

R. S. Bhalerao, J. P. Blaizot, N. Borghini and J. Y. Ol- 

litrault, Phys. Lett. B 627, 49 (2005). 

A. Dumitru, A. Hayashigaki and J. Jalilian-Marian, NucL 

Phys. A 765, 464 (2006); Nucl. Phys. A 770, 57 (2006). 

D. Kharzeev, Y. V. Kovchegov and K. Tuchin, Phys. Rev. 

D 68, 094013 (2003). 

M. A. Braun, Phys. Lett. B 483, 105 (2000). 

R. Baier, A. H. Mueller, D. Schiff and D. T. Son, Phys. 

Lett. B 539, 46 (2002). 

N. Armesto, C. A. Salgado and U. A. Wiedemann, Phys. 

Rev. Lett. 94, 022002 (2005) . 

M. Miller and R. Snellings, nucl-ex/0312008 



S. Manly et al. [PHOBOS ColIaborationJTNucl. Phys. A 
774, 523 (2006)_^ 



S. A. Voloshin, |nucl-th/0606022[ 

R. S. Bhalerao and J. Y. Ollitrault, lnucl-th/0607009', 
C. E. Aguiar, Y. Hama, T. Kodama and T. Osada, Nucl. 
Phys. A 698, 639 (2002). 

Y. Hama, R. P. G. Andrade, F. Grassi, O. So- 
colowski Jr., T. Kodama, B. Tavares and S. S. Padula, 
Nucl. Phys. A 774, 169 (2006); R. Andrade F. Grassi, 
Y. Hama, T. Kodama, O. . J. Socolowski and B. Tavares, 
nucl-th/0511021 

R. Andrade, F. Grassi, Y. Hama, T. Kodama and O. 



Socolowski Jr. nucl-th/0608CI67 
[23] X. 1. Zhu, M. Bleicher and H. Sto ecker, Phys. Rev. C 72, 

064911 (2005); nucl-th/06 0iM9] 
[24] L. V. Gribov, E. M. Levin, and M. G. Ryskin, Pliys. 

Rept. 100, 1 (1983). 
[25] B. B. Back et al. [PHOBOS Collaboration], Phys. Rev. 

C 65, 061901 (2002). 
[26] L. McLerran and R. Venugopalan, Phys. Rev. D 49, 

2233 (1994); ibid. 49, 3352 (1994); 

[27] A. Dumitru and L. D. McLerran, Nucl. Phys. A 700, 492 

(2002) ; A. Dumitru and J. Jalilian-Marian, Phys. Lett. 
B 547, 15 (2002); J. P. Blaizot, F. Gelis and R. Venu- 
gopalan, Nucl. Phys. A 743, 13 (2004). 

[28] A. Krasnitz and R. Venugopalan, Nucl. Phys. B557, 
237 (1999); Phys. Rev. Lett. 86, 1717 (2001); 84, 4309 
(2000); A. Krasnitz, Y. Nara, and R. Venugopalan, Phys. 
Rev. Lett. 87, 192302 (2001); Nucl. Phys. A 727, 427 

(2003) . 

[29] T. Lappi, Phys. Rev. C 67, 054903 (2003). 

[30] A. Krasnitz, Y. Nara and R. Venugopalan, Nucl. Phys. 

A 717, 268 (2003); Phys. Lett. B 554, 21 (2003). 
[31] G. Baym, B. Blattel, L. L. Frankfurt, H. Heiselberg and 

M. Strikman, Phys. Rev. C 52, 1604 (1995). 
[32] D. Molnar and M. Gyulassy, Nucl. Phys. A 697, 495 

(2002), [Erratum-ibid. A 703, 893 (2002)]. 
[33] Z. Xu and C. Greiner, Phys. Rev. C 71, 064901 (2005); 

thep-ph/0509324, 
[34] Y. Nara, Nucl. P hys. A 774, 783 (2006) 
arXiv:nucl-th/05b9052 ; A. Dumitru and Y. Nara, 

Phys. Lett. B 621, 89 (2005) arXiv:hep-ph/0503121]; 

Eur. Phys. J. A 29, 65 (2006) arXiv:hep-ph/0511242'; 

A. Dumitru, Y. Nara and M. Strickland, 

|arXiv:hep-ph/0604149| 



